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Abstract 

Supernovae observations strongly support the presence of a cosmological constant, but its value, 
which we will call apparent, is normally determined assuming that the Universe can be accurately 
described by a homogeneous model. Even in the presence of a cosmological constant we cannot 
exclude nevertheless the presence of a small local inhomogeneity which could affect the apparent value 
of the cosmological constant. Neglecting the presence of the inhomogeneity can in fact introduce 
a systematic misinterpretation of cosmological data, leading to the distinction between an apparent 
and true value of the cosmological constant. We establish the theoretical framework to calculate the 
corrections to the apparent value of the cosmological constant by modeling the local inhomogeneity 
with a ALTB solution. Our assumption to be at the center of a spherically symmetric inhomogeneous 
matter distribution correspond to effectively calculate the monopole contribution of the large scale 
inhomogeneities surrounding us, which we expect to be the dominant one, because of other observations 
supporting a high level of isotropy of the Universe around us. 

By performing a local Taylor expansion we analyze the number of independent degrees of freedom 
which determine the local shape of the inhomogeneity, and consider the issue of central smoothness, 
showing how the same correction can correspond to different inhomogeneity profiles. Contrary to 
previous attempts to fit data using large void models our approach is quite general. The correction 
to the apparent value of the cosmological constant is in fact present for local inhomogeneities of any 
size, and should always be taken appropriately into account both theoretically and observationally. 
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I. INTRODUCTION 

High redshift luminosity distance measurements jl-6] and the WMAP measurements [7, 8| 
of cosmic microwave background (CMB) interpreted in the context of standard FLRW cosmo- 
logical models have strongly disfavored a matter dominated universe, and strongly supported 
a dominant dark energy component, giving rise to a positive cosmological acceleration. 

As an alternative to dark energy, it has been proposed jol, [h]] that we may be at the center 
of an inhomogeneous isotropic universe without cosmological constant described by a Lemaitre- 
Tolman-Bondi (LTB) solution of Einstein's field equations, where spatial averaging over one 
expanding and one contracting region is producing a positive averaged acceleration a^, but 
it has been shown how spatial averaging can give rise to averaged quantities which are not 



observable 



Another more general approach to map luminosity distance as a function of 



redshift Dl(z) to LTB models has been recently proposed 12|, [l3j], showing that an inversion 
method can be applied successfully to reproduce the observed Dl(z). Interesting analysis of 
observational data in inhomogeneous models without dark energy and of other theoretically 



related problems is given for example in 14J-I31] 



Here in this paper we will adopt a different approach. We will consider a Universe with a 
cosmological constant and some local large scale inhomogeneity modeled by a ALTB solution 
[32 1 . For simplicity we will also assume that we are located at its center. In this regard this can 
be considered a first attempt to model local large scale inhomogeneities in the presence of the 
cosmological constant or, more in general, dark energy. Given the spherical symmetry of the 
LTB solution and the assumption to be located at the center our calculation can be interpreted 
as the monopole contribution of the large inhomogeneities which surround us. Since we know 
from other observations such as CMB radiation that the Universe appears to be highly isotropic, 
we can safely assume that the monopole contribution we calculate should also be the dominant 
one, making our results even more relevant. After calculating the null radial geodesies for a 
central observer we then compute the luminosity distance and compare it to that of AC DM 
model, finding the relation between the two different cosmological constants appearing in the 
two models, where we call apparent the one in the ACDM and true the one in ALTB. Our 
calculations show that the corrections to f^ pp , which is the value of the cosmological constant 
obtained from analyzing supernovae data assuming homogeneity, can be important and should 
be taken into account. 



3 



II. LTB SOLUTION WITH A COSMOLOGICAL CONSTANT 



The LTB solution can be written as |33M35| as 

ds 2 = -dt' + + , (1) 

where R is a function of the time coordinate t and the radial coordinate r, E(r) is an arbitrary 
function of r, and f? r = d r R(t,r). The Einstein equations with dust and a cosmological 
constant give 

(R\ 2 _ 2E(r) 2M(r) A 

^yl ~ f? 2 + i? 3 + 3' lJ 

/ \ 2M, r 

p(t ' r)= R^' (3) 
with M(r) being an arbitrary function of r, i? = d t R(t, r) and c = 87rG = 1 is assumed 
throughout the paper. Since Eq. (|2J) contains partial derivatives respect to time only, its 
general solution can be obtained from the FLRW equivalent solution by making every constant 
in the latter one an arbitrary function of r. 

The general analytical solution for a FLRW model with dust and cosmological constant was 



obtained by Edwards [36] in terms of elliptic functions. By an appropriate choice of variables 
and coordinates, we may extend it to the LTB case thanks to the spherical symmetry of both 
LTB and FLRW models, and to the fact that dust follows geodesies without being affected 
by adjacent regions. An anaytical solution can be found by introducing a new coordinate 
rj = rj(t, r) and a variable a by 

I dr)\ r _ 1 
\dt ) r ~R = a 

and new functions by 



(4) 



( ^ 6M ( r ) U \ - 2E ^ (K\ 



Then Eq. (j2J) becomes 



{^) = - Hr)a + — a+ 3 a ' (6) 
where a is now regarded as a function of rj and r, a = a(r),r). It should be noted that the 

coordinate r], which is a generalization of the conformal time in a homogeneous FLRW universe, 

has been only implicitly defined by Eq. (jlj). The actual relation between t and r\ can be obtained 

by integration once a(r], r) is known: 

rv 

t(r),r) = / a(x,r)dx + t b (r) , (7) 
Jo 
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which can be computed analytically, and involve elliptic integrals of the third kind (37)] . 

The function £e(r) plays the role of constant of integration, and is an arbitrary function of r, 
sometime called bang function, since by construction at time t = h(r) we have a(tb(r), r) = 0, 
and correspond to the fact that the big bang initial singularity can happen at different times 
at different positions from the center in a LTB space. In the rest of this paper we will assume 
homogeneous bang, i.e. we will set 

U{r) = 0. (8) 
Inspired by the construction of the solution for the FLRW case get: 

«r,,r) = —, Pj ^K , (9) 

where <f>(x] g 2 , #3) is the Weierstrass elliptic function satisfying the differential equation 
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dx 



- 92(f) - 93 , (10) 



and 



« = Mr) , 92 = \k{rf , 93 = ^ (*k{rf - Ap (r) 2 ) . (11) 

In this paper we will choose the so called FLRW gauge, i.e. the coordinate system in which 
Po(r) is constant. 



III. GEODESIC EQUATIONS AND LUMINOSITY DISTANCE 



We adopt the same method developed in [38j to solve the null geodesic equation written 



in terms of the coordinates (rj,r). Instead of integrating differential equations numerically, we 
perform a local expansion of the solution around z = corresponding to the point (to?0), or 
equivalently (r] , 0), where t = t(r]o, 0). The change of variables from (t, r) to (77, r) permits us 
to have r.h.s. of all equations in a fully analytical form, in contrast to previous considerations 
of this problem which require a numerical calculation of R(t, r) from the Einstein equation (J2j). 
Thus, this formulation is particularly suitable for derivation of analytical results. 

The luminosity distance for a central observer in the LTB space-time as a function of the 
redshift z is expressed as 

D L (z) = (l + z) 2 R(t(z),r(z)) = (l + z) 2 r(z)a(v(z),r(z)) , (12) 
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where \t(z),r(z)j or y(r/(z),r(z)j is the solution of the radial geodesic equation as a function 
of z. The past-directed radial null geodesies is given by 



dt 
dr 



R, r (t,r) 



In terms of z, Eq. ( fl3|) takes the form 391 ] : 



dr 
dz 
dt 

dz 



1 + 2E{r) 



l + 2E(r(z)) 



(13) 



[l + z)R r [r(z),t{z)} ' 
R >r [r(z),t(z)\ 



(14) 



(l + z)R, r [r(z),t(z)] 

The inconvenience of using the (t, r) coordinates is that there is no exact analytical solution 
for R(t,r). So the r.h.s. of Eqs. ( I14jl cannot be evaluated analytically, but we are required 
to find a numerical solution for R first j40j, and then to integrate numerically the differential 
equations, which is quite an inconvenient and cumbersome procedure, and cannot be used to 
derive anaytical results. 



It can be shown 



that in the coordinates (77, r) eqs. (|14p take the form: 
d r t(r),r) + F(rj,r) 



dr) 

dz 
dr 
dz 



(l + z)d v F{rt,r) 
(1 + z)d v F(ri,r) ~ 



p(r),r) 



q(r),r) 



(15) 
(16) 



where 



F(r),r) = 



yjl + 2E(r) yjl - jfc(r)r s 



d r (a(r), r)r) — a d v (a(r], r)r) d r t(r), r) 



(17) 



It is important to observe that the functions p, q, F have explicit analytical forms, making it 
particularly useful do derive anaytical results. 



IV. NUMBER OF INDEPENDENT PARAMETERS AND TAYLOR EXPANSION 
ACCURACY 

In order to find the relation between the apparent and true value of the cosmological constant 
we need in to match the terms in the red-shift expansion : 

D ACDM = jjALTB ^ (lg) 
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Before proceeding in deriving this relation we need to understand clearly how many independent 
parameters we can solve for at different order in the Taylor expansion for Dl(z). After defining 
the expansion of the function k(r) in terms of the dimensionless function K(r): 

k(r) = (a H ) 2 K(r) = K + K x r + K 2 r 2 + .. (19) 

we have 

= j^Ztb , (20) 



Df = m^K^K^.^K^), (21) 

which implies that if we want to match the coefficient A up to order n, we will have a total of 
n + 2 independent parameters to solve for : 

{H^Q^Ko,^,..,^}. (22) 

The matching conditions will imply a constraint over the n + 2 independent parameters, but 
this will not be enough completely determine them, since two of them will always be free. For 
a matter of computational convenience we will choose K , K\ as free parameters and express 
all the other in terms of them. For example from : 



we can get 



from 



we can get 



and in general from 



we can get 



Z?2 — i (23) 

Wa^o^i), (24) 

jjACDM = jjALTB ^ ^ 

W.^ifi), (26) 

jjACDM _ jjALTB ^7) 



^(nr.iffl^i). (28) 



Since our purpose is to find the corrections to the apparent value of the cosmological constant, 
the second order term D 2 is enough. Higher order terms in the redshift expansion will provide 



A(z) 
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FIG. 1: The percentual error A = 100 



,ACDM _j-)ACDM 
Taylor 

D KCDM 



for a third order expansion is plotted as a 



function of the redshift. As it can be seen the error is already quite large at redshift 0.1. Higher order 
expansion does not improve the convergence. 

K2, K%, --Ki-i as functions of {f2 A pp , Kq, K\}, but will not change the analytical relation between 
and Qa which can be derived from eq. ( 123|) . For this reason we will only need the expansion 
up to second order for the luminosity distance. The fact that we have more free parameters 
than constraints implies that the same correction to the apparent value of the cosmological 
constant can correspond to an infinite number of different inhomogeneity profiles. 

The corrections we calculate are accurate within the limits of validity of the Taylor expansion 
^Taylor ■ ^ turns out that in the flat case we consider the error is quite large already at a redshift 
of about 0.2 as shown in the figure. This implies that the corrections should also be valid only 
within this low redshift range, since even if we are exactly matching the coefficients, the Taylor 
expansion of the AC DM best fit formula itself is not very accurate. This could be overcome 
by implementing other types of expansions or numerical methods, such as Pade for example, 
with better convergence behavior, but we'll leave this to a future work. 

V. CENTRAL BEHAVIOR 

A function of the radial coordinate f(r) is smooth at the center r = only if all its odd 
derivatives vanish there. This can be shown easily by looking at the partial derivatives of even 
order of this type for example: 




(29) 



S 



where {x, y, z] are the cartesian coordinates related to r by r 2 = x 2 + y 2 + z 2 . Quantities of 
the type above diverge at the center if d 2m+1 f(r) ^ for 2m + 1 < In. If for example the first 
derivative f'(0) is not zero, then the laplacian will diverge. This implies that including linear 
terms expansions for k(r) and t h {r) we are considering models which are not smooth at the 
center. The general central smoothness conditions are: 

k 2m+ i = 0, (30) 
tl m+1 = 0, (31) 
2m + 1 < i, (32) 

which must be satisfied for all the relevant odd powers coefficients of the central Taylor expan- 
sion. In our case this implies that if we only want to consider centrally smooth inhomogeneities 
then we need to set to zero all the odd derivatives of K(r) 

K 2m+1 = (33) 

The consequence of this smoothness conditions is that the exact matching of the Taylor expan- 
sion is possible only up to order five when we have five constraints equations 

D ACDM = D ALTB 1<*< 5 , (34) 

and five free parameters 

H£ LTB ,n A ,K ,K 2 ,K 4 (35) 

implying there is a unique solution. Going to higher order there will be more equations than 
free parameters making the inversion problem impossible. This means that the effects of a 
different value of the cosmological constant cannot be mimicked by a smooth inhomogeneity, 
as far as the exact matching of the Taylor expansion is concerned. From a data analysis point 
of view this limitation could be easily circumvented, since these considerations are based on 
matching the Taylor expansion of the best AC DM fit, which is quite different from fitting 
the actual data. Also it turns out that the Taylor expansion Dj,%$Jf(z) is more accurate at 
second order than at any other order as shown in the figure, implying that exact matching 
beyond second order is practically irrelevant from a data fitting point of view. Under these 
considerations the inversion problem can be considered still effectively undetermined since by 
matching up to second order we have two equations and three parameters: 

H ALTB ,Q A ,K (36) 



For completeness of the analysis we mention that after counting the number of independent 
parameters we can easily conclude that the inversion problem remain undetermined for the 
third order, and has a unique solution for the fourth and fifth order as shown above. 



VI. CALCULATING THE LUMINOSITY DISTANCE 

In order to obtain the redshift expansion of the luminosity distance we need to use the 
following: 

k(r) = (a H ) 2 K(r) = K + K x r + K 2 r 2 + .. (37) 

t(rj,r) = bo(rj) + h {rj)r + b 2 (r])r 2 + •• (38) 

It should be noted that linear terms will in fact lead to central divergences of the laplacian in 
spherical coordinates, which correspond to a central spike of the energy distribution 29_|, [30] , 
but an appropriate local averaging of the solution can easily heal this behavior, and we include 
them here because they give the leading order contribution. Since we are interested in the 
effects due to the inhomogeneities we will neglect k in the rest of the calculation because this 
corresponds to the homogeneous component of the curvature function k(r). 



Following the same approach given in 



we can find a local Taylor expansion in red-shift 



for the geodesies equations, and then calculate the luminosity distance: 



D 



ALTB, 



n ALTB _ - 

1 ~ iV 



:i + z) 2 r(z)a ALTB (r)(z),r(z)) = D* LTB z + Dp™-- 
1 



f +D ALTBJi 



Z° + 



(39) 



D 



ALTB 



545 1 (fi^ ue - l) 2 + 18^(0^ - 1) - 18/i , r ( fi 



^true^2 



+30/i ,r^T e - 12 V + GK^X^ - lOifi + 27(ftf e ) 2 



where we have introduced the dimensionless quantities K , Ki, B\, B[, /i , r according to 



/ d t ,a(t,r) 
Bi(v) = Hv)^ 1 , 



t=t ,r=0 



' d v a(r),r) 
, a(r?,r) 2 



77=770 ,r=0 



B 1 

ho,r 

to 



bi(vo)a \ 
dB x {ri) 



drj 
1 



(a i?o 



v=vo 
d r a(r), r) 



a H a(r],r) 



77=770, 7*=0 



(40) 

(41) 

(42) 
(43) 
(44) 

(45) 
(46) 
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and used the Einstein equation at the center (77 = t/q, r = 0) 

1 = 14(0) + tt M + = -K + tt M + Oa, (47) 

= (48) 

n « = 4?' (49) 

f!/v = 4- (50) 

Because of our coordinate choice O w is independent of r, and all the radial dependence goes 
into fifc(r). Note that apart from the central curvature term K , the inhomogeneity of the LTB 
space is expressed in ho >r , which encodes the radial dependence of the scale factor. Details of 
these rather cumbersome calculations are provided in a separate companion paper, but it should 
be emphasized that in order to put the formula for the luminosity distance in this form it is 
necessary to manipulate appropriately the elliptic functions and then re-express everything in 
terms of physically meaningful quantities such as H . 

VII. CALCULATING D L (z) FOR ACDM MODELS. 

The metric of a ACDM model is the FLRW metric, a special case of LTB solution, where : 

Po(r) oc const, (51) 

k(r) = 0, (52) 

t b (r) = 0, (53) 

a(t,r) = a(t). (54) 

We will calculate independently the expansion of the luminosity distance and the redshift 
spherical shell mass for the case of a flat ACDM, to clearly show the meaning of our notation, 
and in particular the distinction between Vt a £ p and fi^" 6 . We can also use these formulas to 
check the results derive before, since in absence of inhomogeneities they should coincide. 
One of the Einstein equation can be expressed as: 



h acdm (z) = # y (1 - (^) + = H oy /(i-nT')(i + z) 3 + nf p . (55) 

We can then calculate the luminosity distance using the following relation, which is only valid 
assuming flatness: 

r)KCDM ( \ /-1 1 ^\ [ Z dz' t^ACDM , nACM 2 , nACDM 3 _j_ 

L {Z) = i 1 + Z) J o JJACDM^ =D l Z + D 2 Z + D 3 Z + - ( 56 ) 
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From which we can get: 



D ACDM = 1 (5?) 
n ACDM _ 3^/ P + 1 

^ 2 " 4# ' ( j 

We can check the consistency between these formulae and the ones derived in the case of LTB 

by setting: 

Kt = B 1 = B[ = K = V = , (59) 

(60) 

which corresponds to the case in which = f^ ue . 

VIII. RELATION BETWEEN APPARENT AND TRUE VALUE OF THE COSMO- 
LOGICAL CONSTANT 

So far we have calculated the first two terms of the redshift expansion of the luminosity 
distance for ALTB and AC DM model. Since we now that the latter provides a good fitting for 
supernovae observations, we can now look for the ALTB models which give the same theoretical 
prediction. From the above relations we can derive : 

H ALTB = H£ CDM , (61) 
T = 2 7 (^1 _ 1} 54^(1^-1085!^ + 5AB 1 + 18B[n t r-18B' 1 
-18/i ,r(^A Me ) 2 + 30V^A Ue - 12 V + QK^X^ - lOi^i 
+27tt% ue {tt% ue - 1) 



1 



(62) 



nX ue 



1 



6(6Bi - 2h , r + 3) 



^(36£i - 6B[ - 10h 0:r - 2K X + 9nf p + 9) 2 + 



\ 1/2 

-4(6^ - 2h 0>r + 3)(54 J B 1 - 185; - 12V ~ 10K i + 27^ pp ) j - 36^ 
+QB[ + 10V + 2K i ~ 9 ( fi A P - !) 



(63) 



We can also expand the above exact relations assuming that all the inhomogeneities, can be 
treated perturbatively respect to the AC DM , i.e. {Ki, Bi, B[} oc e, where e stands for a small 
deviation from FLRW solution : 



Q true = n ap P _ ^ ^ (21 B^T ~ VT + ^("T " 1) " 9V("T) + ^V^T 

-6V + ^K 1 Q^ P - hK x ) + 0(e 2 ) . (64) 
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As expected all these relations reduce to 

in the limit in which there is no inhomogeneity, i.e. when K\ — B\ — B[ — ho jr = 0. 

IX. CONCLUSIONS 

We have derived for the first time the correction due to local large scale inhomogeneities 
to the value of the apparent cosmological constant inferred from low redshift supernovae ob- 
servations. This analytical calculation shows how the presence of a local inhomogeneity can 
affect the estimation of the value of cosmological parameters, such as This effects should 
be properly taken into account both theoretically and observationally. By performing a local 
Taylor expansion we analyzed the number of independent degrees of freedom which determine 
the local shape of the inhomogeneity, and consider the issue of central smoothness, showing how 
the same correction can correspond to different inhomogeneity profiles. We will address in a 
future work the estimation of the magnitude of this effect based on experimental bounds which 
can be set on the size and shape of a local inhomogeneity and the fitting of actual supernovae 
data. It is important to underline here that we do not need a large void as normally assumed in 
previous studies of LTB models in a cosmological context. Even a small inhomogeneity could 
in fact be important. 

In the future it will also be interesting to extend the same analysis to other observables 
such as barionic acoustic oscillations (BAO) or the cosmic microwave background radiation 
(CMBR), and we will report about this in separate papers. Another direction in which the 
present work could be extended is modeling the local inhomogeneity in a more general way, for 
example considering not spherically symmetric solutions. From this point of view our calcula- 
tion could be considered the monopole contribution to the general effect due to a local large 
scale inhomogeneity of arbitrary shape. Given the high level of isotropy of the Universe shown 
by other observations such as the CMB radiation, we can expect the monopole contribution we 
calculated to be the dominant one. 

While this should be considered only as the first step towards a full inclusion of the effects 
of large scale inhomogeneities in the interpretation of cosmological observations, it is important 
to emphasize that we have introduced a general definition of the concept of apparent and true 
value of cosmological parameters, and shown the general theoretical approach to calculate the 
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corrections to the apparent values obtained under the standard assumption of homogeneity. 
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